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Abstract 

The H-coalescent processes were initially studied by Mohle and Sagitov 
(2001), and introduced by Schweinsberg (2000) in their full generality. They 
arise in the mathematical population genetics as the complete class of scal- 
ing limits for genealogies of Cannings' models. The E-coalescents generalize 
A-coalescents, where now simultaneous multiple collisions of blocks are pos- 
sible. The standard version starts with infinitely many blocks at time 0, 
and it is said to come down from infinity if its number of blocks becomes 
immediately finite, almost surely. This work builds on the technique intro- 
duced recently by Berestycki, Berestycki and Limic (2009), and exhibits a 
deterministic "speed" function - an almost sure small time asymptotic to 
the number of blocks process, for a large class of H-coalescents that come 
down from infinity. 
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1 Introduction 



Kingman's coalescent [T5| [16] is one of the central models of mathematical popu- 
lation genetics. From the theoretical perspective, its importance is linked to the 
duality with the Fisher- Wright diffusion (and more generally with the Fleming- 
Viot process). Therefore the Kingman coalescent emerges in the scaling limit of 
genealogies of all evolutionary models that are asymptotically linked to Fisher- 
Wright diffusions. From the practical perspective, its elementary nature allows 
for exact computations and fast simulation, making it amenable to statistical 
analysis. 

Assume that the original sample has m individuals, labeled {1, 2, . . . , m}. One 
can identify each of the active ancestral lineages, at any particular time, with a 
unique equivalence class of {1,2,..., m} that consists of all the individuals that 
descend from this lineage. In this way, the coalescent event of two ancestral 
lineages can be perceived as the merging event of two equivalence classes. Ignoring 
the partition structure information, one can now view the coalescent as a block 
(rather than equivalence class) merging process. 

Kingman's coalescent corresponds to the dynamics where each pair of blocks 
coalesces at rate 1. Hence, if there are n blocks present in the current configura- 
tion, the total number of blocks decreases by 1 at rate (™) • Using this observation 
and elementary properties of exponential random variables, one can quickly con- 
struct the standard version of the process, which "starts" from a configuration 
containing an infinite number of individuals (particles, or blocks) at time 0, and 
has the property that its configuration consists of finitely many blocks at any 
positive time. 

The fact that in the Kingman coalescent dynamics only pairs of blocks can 
merge at any given time makes it less suitable to model evolutions of marine 
populations or viral populations under strong selection. In fact, it is believed 
(and argued to be observed in experiments, see e.g. [IT]) that in such settings 
the reproduction mechanism allows for a proportion of the population to have the 
same parent (i.e., first generation ancestor). This translates to having multiple 
collisions of the ancestral lineages in the corresponding coalescent mechanism. 

A family of mathematical models with the above property was independently 
introduced and studied by Pitman [21] and Sagitov [22] under the name A- 
coalescents or coalescents with multiple collisions. Almost immediately emerged 
an even more general class of models, named H- coalescents or coalescents with 
simultaneous multiple collisions or exchangeable coalescents. The Greek letter S 
in the name is a reference to the driving measure H (see Sections 12.21 and 12.41 for 
details). The S-coalescent processes were initially studied by Mohle and Sagitov 
[T9] , and introduced by Schweinsberg [21] in their full generality. In particular, it 
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is shown in [19] that any limit of genealogies arising from a population genetics 
model with exchangeable reproduction mechanism must be a S-coalescent. 



The current paper uses the setting and several of the results from [24J that will 

be recalled soon. Formally, under the S-coalescent dynamics, several families of 
blocks (with two or more blocks in each family) may (and typically do) coalesce 
simultaneously. The S-coalescents will be rigorously defined in the next section. 

More recently, Bertoin and Le Gall [9j established a one-to-one correspondence 
between a class of processes called stochastic flows of bridges and S-coalescents, 
and then constructed in [10] the generalized Fleming- Viot (or A-Fleming-Viot) 
processes that have A-coalescent processes as duals; Birkner et al. p2] recently ex- 
tended this further by constructing for each driving measure S the "E-Fleming- Viot 
processes, dual to the corresponding S-coalescent process; Durrett and Schweins- 
berg [TJ] showed that genealogies during selective sweeps are well-approximated 
by certain S-coalescents; and Birkner and Blath [TT] initiated a statistical study 
of coalescents with multiple collisions. 

Generalizations of A-coalescents to spatial (not a mean-field) setting are stud- 
ied by Limic and Sturm [18], and more recently by Angel et al. [T] and Barton et 
al. [2] . The reader can find detailed information about these and related research 
areas in recent texts by Berestycki [3] and Bertoin [8]. 

Let N~ = N := (N(t), t > 0) be the number of blocks process correspond- 
ing to a particular standard (meaning lim t ^ + N(t) = oo) S-coalescent process. 
Moreover, suppose that this S-coalescent comes down from infinity, or equiva- 
lently, assume that P(N(t) < oo, Vt > 0) = 1 (see the end of Section I2T41 for a 
formal discussion). From the practical perspective, it seems important to under- 
stand the nature of the divergence of N(t) as t decreases to (see [5] for further 
discussion and applications). 

The main goal of this work is to exhibit a function i>n = v : (0, oo) — > (0, oo) 
such that 



We call any such v the speed of coming down from infinity (speed of CDI) for 
the corresponding S-coalescent. The exact form of the function v is implicit and 
somewhat technical, see Theorems [1] or [10] for a precise statement. Moreover, the 
speed is obtained under (relatively weak) additional "regularity" condition. 

The coming down from infinity property was already studied by Schweinsberg 
[24] in detail. The speed of coming down from infinity for general S-coalescents 
has not been previously studied. In Berestycki et al. [3] the speed of CDI of any 
A-coalescent that comes down from infinity was found using a martingale-based 
technique. A modification of this technique will be used presently to determine 





1, almost surely. 



(1) 
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the above v, and the steps in the argument that carry over directly to the current 
setting will only be sketched. 

In the A-coalescent setting, weaker asymptotic results (than (CQ)) on N~fv = 
N A /v can be deduced by an entirely different approach, based on the theory of 
Levy processes and superprocesses. This link was initially discovered in [HI [7] 
in the special case of so-called Beta-coalescents, and recently understood in the 
context of general A-coalescents in (5]. It is worthwhile pointing out, that for any 
"true" H-coalescent (meaning that simultaneous multiple collisions are possible in 
the dynamics), an approach analogous to [5j seems rather difficult to implement 
(to start with, the expression ( 120|) . unlike ( 1T5]) . does not seem to be directly linked 
with any well-known stochastic process). Indeed, the martingale technique from 
[I] has at least three advantages: (i) it yields stronger forms of convergence, more 
precisely, ([1]) and its counterparts in the L p -sense, for any p > 1 (cf. |4J Theorems 1 
and 2); it yields explicit error estimates needed in the frequency spectrum analysis 
(see jH E] for details); and (iii) it extends to the S-coalescent setting as will be 
explained shortly. 

It is not surprising that the form of the "candidate" speed of CDI for the 
H-coalescent is completely analogous to that for the A-coalescent. What may be 
surprising is that there are S-coalescents that come down from infinity but their 
candidate speed is identically infinite. And also that there might be coalescents 
that come down from infinity but faster than their finite candidate speed. (Remark 
|6] in Section 12.41 explains how neither of these can occur under a A-coalescent 
mechanism.) The question of whether an asymptotic speed still exists in such 
cases remains open. This discussion will be continued in Section 13.21 

The rest of the paper is organized as follows: Section [2] introduces various 
processes of interest (including in Section 12.61 the novel "color-reduction" and 
"color-joining" constructions that might be of independent interest) and presents 
a "preview" of the main result as Theorem [TJ Section [3] contains a "matured" 
statement of the main result (Theorem [TO]) . followed by a discussion of some 
of its immediate consequences, and of the significance of a certain "regularity 
hypothesis" , while Section @] is devoted to the proof of Theorem [TUJ 

2 Definitions and preliminaries 
2.1 Notation 

In this section we recall some standard notation, as well as gather less standard 
notation that will be frequently used. 

Denote the set of real numbers by K. and set M + = (0, oo). For a, 6 6 R, denote 
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by a A h (resp. a V h) the minimum (resp. maximum) of the two numbers. Let 

A := {{x 1 ,x 2 , . . .) : x x > x 2 > . . . > 0, < 1}, (2) 

i 

be the infinite unit simplex. If x = (x x , x 2 , . . .) G A and c G M, let 

CX = (cXi, CX 2 , • • •)• 

Denote by the zero (0,0...,) in A. 

Let N := {1,2, . . .}, and V be the set of partitions of N. Furthermore, for 
n eN denote by V n the set of partitions of [n] := {1, . . . , n}. 

If / is a function, defined in a left-neighborhood (s — e, s) of a point s, denote 
by f(s—) the left limit of / at s. Given two functions /, g : K + — > R + , write 
/ = 0(g) if limsup f(x)/g(x) < oo, / = o(g) if limsup f(x)/g(x) = 0, and 
/ ~ g if lim/(x)/^(x) = 1. Furthermore, write / = Q(g) if both / = 0(g) and 
g = 0(f). The point at which the limits are taken is determined from the context. 

If J 7 = (Ft, t > 0) is a filtration, and T is a stopping time relative to J 7 , denote 
by Tt the standard filtration generated by T, see for example [13], page 389. 

For v a finite or a- finite measure on A or on [0,1], denote the support of v by 
supp(v). 



2.2 S-coalescents 

Let H be a finite measure on A, and write 

H = S + a<5o, 

where a > and S ((0,0, . . .)) = 0. As noted in [21], we may assume without 
loss of generality that S is a probability measure. The S-coalescent driven by the 
above S is a Markov process (II t ,i > 0) with values in V (the set of partitions of 
N), characterized in the following way. If n G N, then the restriction (Ilf 1 , t > 0) 
of (Ht,t > 0) to [n] is a Markov chain, taking values in V n , such that while Il[ n ^ 
consists of b blocks, any given fci-tuple, k 2 -tup\e,. . ., and fc r -tuple of its blocks 
(here YH=i ki <b and fcj > 2, i — 1, . . . , r) merge simultaneously (each forming 
one new block) at rate 



A^fci,...,^;! — / 2 C^QXJ, 



where s := 6 — X^ =1 fci is the number of blocks that do not participate in the 
merger event, and where the sum Y\ in the above summation stands for 
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the infinite sum V°° , Y]°° , . , . ... Y]°° ^ . , over r + I different 

indices. It is easy to verify that each such coalescent process has the same rate of 
pairwise merging 

A 2;2;0 = S(A) = 1. (3) 



2.3 Preview of the small-time asymptotics 

One can now state the central result of this paper. Given a probability measure 
5 as above, for each t > denote by N~(t) the number of blocks at time t in the 
corresponding (standard) S-coalescent process. Define 

Ml) ■= I ^^i' 1 ^^ 5(dk), q > 0, 
J A l^i=l x i 



and 



v s (t) := inf Is > : f dq < t\ , t > 0. 

{ Js ^a(g) J 



Theorem 1. // &o£A 

S({x G A : Xi = 1 for some finite n}) = 



i=i 



and 



/ ^ S(&) < oo, 

J A 2-a=l X i 



lim — -r-r- = 1, almost surely, 
t^o+ v s (t) 

where oo/oo = 1. In particular, under the above assumptions, the quantity Vs(t) 
is finite for (one and then for) all t > 0, if and only if the H- coalescent comes 
down from infinity. 

Most of the sequel is devoted to explaining the above implicit definition of the 
speed Vs, as well as the significance of the two hypotheses in Theorem [TJ 

As already mentioned, Theorem [I] is restated as Theorem [10] in Section [3j 
which is proved in Section HI The additional condition S({0}) = in Theorem [TD] 
is not really restrictive, since the case where S({0}) > is already well-understood 
(cf. Remark [3] below). 
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2.4 Basic properties of S-coalescents 



Recall the setting and notation of Section 12.21 

If supp(H) C {(x, 0, 0, . . .) : x G [0,1]}, the resulting S-coalescent is usually 
called the A-coalescent, where A is specified by 

A{dx) :=E{d{x,0,...)). (4) 

The transition mechanism simplifies as follows: whenever 11^ consists of b blocks, 
the rate at which any given fc-tuple of its blocks merges into a single block equals 

\ b ,k= [ x k - 2 (l-x) b - h A(dx). (5) 

J '[0,1] 

Note that mergers of several blocks into one are still possible here, but multiple 
mergers cannot occur simultaneously. 

We recall several properties of the H-coalescents carefully established in |24j . 
the reader is referred to this article for details. The H-coalescents can be con- 
structed via a Poisson point process in the following way. Assume that 5 = S , 
or equivalently that S does not have an atom at (see also Remark [3] below). Let 

*(■) = £ WO (6) 

ken 

be a Poisson point process on R + x A with intensity measure dtCg)H(dx)/ x l- 
Each atom (t, x) of 7r influences the evolution of the process IT as follows: to each 
block of II(t— ) assign a random "color" in an i.i.d. fashion (also independently of 
the past) where the colors take values in NU (0, 1) and their common distribution 
P x is specified by 

oo 

P x ({i}) =x h i>l and P x (du) = (1 - a*) du, u e (0, 1); (7) 

i=i 

given the colors, merge immediately and simultaneously all the blocks of equal 
color into a single block (note that this can happen only for integral colors), while 
leaving the blocks of unique color unchanged. 

Note that in order to make this construction rigorous, one should first consider 
the restrictions (U^(t),t > 0), since the measure S(c?x)/ x 1 ma y have (and 
typically will have in the cases of interest) infinite total mass. Given a fixed time 
s > 0, a small s > and any n £ N, it is straightforward to run the above 
procedure using only the finite number of atoms of 7r that are contained in [0, s] x 
{x G A : Y^iLi x i > £ }- Denote the resulting process by (IT n )' e (t), t G [0, s}). The 



7 



following observation is essential: the atoms of it contained in the "complement" 
[0, s] x {x G A : — £ }i together with the coloring procedure, influence 

further the state of (Jl^' £ (t),t 6 [0, s]) on the event A e]S , n (by causing additional 
mergers during [0,s]), where 

PU„„) < 1 - exp ( fdt / _^!^S(dx) 



. i=l X i 

and where the right-hand-side goes to as £ — >■ 0. Indeed, i? x (n) is the probability 
that under P x (assuming n blocks present at time t—) at least two of the blocks 
are colored by the same color. For the benefit of the reader we include the exact 
expression for this probability: 

oo oo 

n-1 



Rx(n) = 1 - (1 - Xj) n - n Xj(l - Xj 

i=l i=l 

(\ oo oo 
2 ) ^ ] ^ ] x ii x ii.i}- — x ii ~ x 
ii = l i 2 = l,i2=£ii 



00 00 n— 1 n— 1 

11=1 in-l=l,in-10{H,i2,— ,i»-2} £=1 €=1 

oo oo n 

- E- £ II 

*1=1 J„ = l,l„^{ll,J2v>«n-l} 

Note that R^{n) < (™) so that 

tZt / to^2 S ( rfx ) ^ s ( 1 / ~(dx) ^ 0, as £ ^ 0. 



, i=i a: < W Mn{E£^<£} 

In this way one obtains a coupling (that is, a simultaneous construction on a single 
probability space) of the family of processes (n( n )- e (t), t e [o,s]), as £ e (0,1), 
and can define as the limit IT( n ) := lim E _K) n' n ) ,£ on [0,s]. Moreover, the above 
construction is amenable to appending particles/blocks to the initial configuration, 
hence it yields a coupling of 

(fl (n) ' £ (t), t e [0, s]), as £ e (0, 1) and fiGN. (8) 

An interested reader is invited to check (or see [21]) that the limit 

IT := lim IT (n) (9) 
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is a well-defined realization of the S-coalescent, corresponding to the measure 5. 
We will denote its law simply by P (rather than by Pg). 

Consider the above S-coalescent process II. Let E = {N(t) = oo for all t > 0}, 
and F = {N(t) < oo for all t > 0}. 

Definition 2. We say that a H-coalescent comes down from infinity if P(F) = 1. 
Let 

n 

Af : = {x G A : 2_^ Xi = 1 for some finite n}. (10) 

i=i 

Lemma 31 [21] extends Proposition 23 of Pitman [21J to the S-coalescent setting. 
It says that provided S(Aj) = 0, there are two possibilities for the evolution of 
N: either P(E) = 1 or P(F) = 1. 

Remark 3. A careful reader will note that the above Poisson point process (PPP) 
construction assumed 5({0}) = 0. It is possible to enrich it with extra pairwise 
mergers if S({0}) > 0, see [21] for details. For the purposes of the current study 
this does not seem to be necessary. Indeed, by the argument of [1] Section 4.2, 
one can easily see that if S((0, 0, . . . , )) = a > 0, then the corresponding 5- 
coalescent comes down from infinity, and moreover its speed of CDI is determined 
by a. More precisely, such a S-coalescent comes down faster than the A-coalescent 
(II a (s),s > 0) corresponding to A(dx) = aS (dx) (note that Il a is just a time- 
changed Kingman coalescent), and slower than (n o ((l +e)s),s > 0), for any 
e > 0. 

In the rest of this paper we will assume that S({0}) = 0, or equivalently, that 

Remark 4. The condition S(A/) = is similar (but not completely analogous, 
see next remark) to the condition A({1}) = for A-coalescents. It is not difficult 
to construct a H-coalescent, such that S(A/) > and P{E) = P(F) = 0. Take 
some probability measure 5' on A such that the corresponding S-coalescent does 
not come from infinity, and define 5 = (1 — a)S' + au, for some a G (0, 1) and 
some probability measure v on Af (for example v(dx) = 5(i/2,i/2,o,...)(^ x ))- Then 
its block counting process stays infinite for all times strictly smaller than T*, and 
it is finite for all times larger than or equal to T*, where 

T* := inf{s : tt({s} x A/) > 0} 

has exponential (rate a) distribution, hence is strictly positive with probability 1. 
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Remark 5. It may be surprising that there are measures S satisfying S(A/) = 1, 
and such that the corresponding S-coalescent comes down from infinity. Note 
that there is no analogy in the setting of A-coalescents, since if A({1}) = 1 (and 
therefore A([0, 1)) = 0), the only such "A-coalescent" will contain a single block 
for all times. 

It was already observed by Schweinsberg j2l] Section 5.5 that if the quantity 




is infinite, the corresponding S-coalescent comes down from infinity. Moreover, 
if the above quantity is positive and finite, the corresponding S-coalescent comes 
down from infinity if and only if the S-coalescent corresponding to S'(c/x) = 
H(<ix)l{ xe A\A / } comes down from infinity, and the speed of CDI is determined 
by S'. This type of coalescent was already mentioned in Remark HI The reader 
should note that if such a S-coalescent does not come down from infinity, then 
P(E) = P(F) = 0. 

Henceforth we will mostly assume that S(A/) = 0. 



2.5 Coming down from infinity revisited 

In this section we assume that S(Ay) = 0, as well as S({0}) = 0. A sufficient 
condition for a S-coalescent to come down from infinity was given by Schweinsberg 
[24] . For ki, i = 1, . . . ,r such that ki > 2 and s := b — YH=i k — 0> define 
N(b; ki, . . . , k r ; s) to be the number of different simultaneous choices of a fci-tuple, 
a A;2-tuple,. . . and a fc r -tuple from a set of b elements. The exact expression for 
N(b; ki, . . . , k r ; s) is not difficult to find (also given in [24J display (3)), but is not 
important for the rest of the current analysis. Let 

[b/2} 

lb ■= ^2 ^ ( b - r - s)N(b;h, . . . ,k r ;s)\ b]ku ..., kr] s 

r=l {ki,...,k r } 

be the total rate of decrease in the number of blocks for the S-coalescent, when 
the current configuration has precisely b blocks. 

Given a configuration consisting of b blocks and an x G A, consider the coloring 
procedure (J7|), and define 

b 

y( b ) l{ith block has color £}, ^ £ N, (11) 

J'=l 
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so that has Binomial(6, xi) distribution. Due to the PPP construction of the 
previous subsection, we then have 

,(&) 



7b = / ^oo 2 "(" X 



A 



^'r-" ^ (12) 

Proposition 32 in [21] says that the S-coalescent comes down from infinity if 

oo 
6=2 

Let £ G (0, 1) be fixed. Recall (j2J) and define 

A e := {x G A : < 1 -e}. (14) 

Proposition 33 in [21] says that (TIB"]) is necessary for coming down from infinity if 
also 

r l 

/ — - H(dx) < oo, for some e > 0. 

Moreover [21j provides an example of a S-coalescent that comes down from infinity, 
but does not satisfy ( fl~3l) . More details are given in Section I3~2l 

The CDI property for A-coalescents is, in comparison, completely understood. 
Define 

MQ) ■= I {e ~ qX ~] + qX) Kdx), (15) 

J [0,1] x 

and note that 7& simplifies to Ylk=2(^ ~ -0 (fc)-V fc > with as in ([SJ. The original 
sharp criteria is due to Schweinsberg [23]: a particular A-coalescent comes down 
from infinity if and only if (Tl3|) holds. Bertoin and Le Gall [10] observed that 

7b = e(Mb)), (16) 

and that therefore the CDI happens if and only if 

< oo, (17) 



for some (and then automatically for all) a > 0. 
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Remark 6. A variation of the argument from Berestycki et al. JTJ provides 
an independent (probabilistic) derivation of the last claim. More precisely, let 
N n (t) = #II (ri) (t), t > and let v n be the unique solution of the following Cauchy 
problem 

v '{t) = -ip(v(t)), v(0) = n. 

Use the argument of jl] Theorem 1 (or see Part I in Section 14.21 for analogous 
argument in the S-coalescent setting) to find n < oo, a G (0, 1/2), and C < oo 
such that 

n>„ V.-™ ""W 

happens with overwhelming (positive would suffice) probability, uniformly in small 
s. Finally, note that v 11 satisfies the identity 

dq s, s > 0, (18) 



n 



sup 

te[o,s] 



therefore \im n v n {s) < oo if and only if ( ITT)) holds. 

Moreover, it was shown in jl] that under condition (fTTj) . a speed 1 1— y v(t) of CDI 
is specified by 

f4 = U>0- (19) 

Jv(t) WKW 



Consider again the general S-coalescent setting. In analogy to (fT5l) . define 



m = mq) ■= I ^ =ll Voo 2 ; ( 2 o) 



Note that the above integral converges since 

e- z - 1 + z < z 2 /2, for all z > 0, (21) 

and so in particular ipziq) < <? 2 /2, for any probability measure H on A. It is easy 
to check that q V ; h( ( ?) is an infinitely differentiable, strictly increasing, and 
convex function on R + , as well as that ipsiq) ~ ? 2 /2 as q — > 0. Therefore, if 

rfg <oo (22) 



for some a > 0, the same will be true for all a > 0, and irrespectively of that 
Jo dq/ip{q) = oo, for any a > 0. 
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Lemma 7. The function q h-> ip(q)/q is strictly increasing. 

The proof (straightforward and left to the reader) is analogous to that for jl] 
Lemma 9. 

Lemma 8. The conditions ( TJ3j) and |Hj] are equivalent. 

Proof. It suffices to show the order of magnitude equivalence ffTBl) in the current 
setting. Use expression ( fl2l) for 75. Note that if x G [0, 1] and b > 1, then 
e -&* > (1 - in fact for x G [0, 1) 

e~ bx - (1 - a;) 6 = e~ bx (l - exp j-6 f y + y + . 

If a; > 1/4 and 6 > 16 it is clearly true that e~ bx — (1 — x) b < bx 2 . For a; < 1/4 
we have Ylj>2 x ^ — 4x 2 /3, and so 

1 - exp I -6 f y + y + . . )j I < 1 - exp (-jj^ 2 ) < hx 2 . 

We conclude that 

< e ~ bx - (1 - x) b < bx 2 , for all 6 > 16 and x G [0, 1]. (23) 
As a consequence, 



DC 



f=l ^=1 i=l 

where 0(Y^i G [— x « 2 ' 0]- By integrating over S(dx)/ x i> we § et 
7fe = ^(6) + 0(6), for some 0(6) G [-6, 0], 

implying 7 b = 0(^(6)). It is easy to check directly from (fT2"j) that 7f, + i — 7& > 
7& — 7&-i > 0, for any 6 > 3 implying 6 = 0(t&). Using convexity of ip, we now 
have that either ip{b) = 0(6), or 6 = o(ip(b)) so that 

''' l + o(l). 



^(6) 

In both cases we have 'jb = ®(ip(b)). □ 
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Assuming (122]) (or equivalently, (TT3T) ) . one can define 



u s (t) 




dq 



G R+, t>0 



(25) 



and its cadlag inverse 



OO 



1 



Us(i) = v(t) := inf <^ s > : 



dq <t> , t > 0. 



(26) 



Call thus defined the candidate speed. In fact, due to the continuity and strict 
monotonicity of u, Un is again specified by (1T9]) . with replacing If ([!]) holds 
with f = f=, we will sometimes refer to the candidate speed i>« as the irue speed 



Note that (1261) makes sense regardless of (1221) . and yields t>~(£) = oo, for each 
£ > 0, if (and only if) (122]) fails. We will say that that the S-coalescent "has an 
infinite candidate speed" in this setting. 

Due to the fact ips(q) < <? 2 /2 (cf. discussion following ( 120]) ) we have 

Corollary 9. If (TJj holds, then v E (t) > 2/t, for t > 0. 

One could try to rephrase the corollary by saying that among all the S- 
coalescents (satisfying ([3])), the Kingman coalescent is the fastest to come down 
from infinity at speed t \-> 2/t (as is the case in the A-coalescent setting, cf. jl] 
Corollary 3). However, there are examples of H-coalescents with infinite candidate 
speed that do come down from infinity. Moreover, there are coalescents that come 
down from infinity, and that have finite candidate speed i>s, but the methods of 
this article break in the attempt of associating N = iV = and Vs at small times. 
The existence of a (deterministic) speed, and its relation to the function t > 2/t 
in these situations, are open problems. Vaguely speaking, such "difficult cases" 
correspond to measures H for which there exists a set A*? ~ A/ such that 



For rigorous statements see Section [3j 

2.6 Two operations on S-coalescents 

In this section we consider two variations of the PPP construction ([H])-©, each 
of which gives a probabilistic coupling of the original S-coalescent with a simpler 
S-coalescent. 



of GDI. 
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Given a realization of the Poisson point process ([6]) and the coloring of (J7J), 
define the 5-reduction or (5 -color-reduction) U r s (r stands for "reduction") of II 
to be the partition valued process constructed as follows: immediately after each 
coloring step (and before the merging) run Bernoulli (5) random variable for each 
block, independently over the blocks and the rest of the randomness, and for each 
of the blocks having this new value 1, resample its color from the uniform U[0, 1] 
distribution, again independently from everything else. Note that the "reduction" 
in the name refers to reducing the coloring (atom) weights, however this has the 
opposite effect on the number of blocks. Indeed, the above procedure makes some 
blocks that share (integer) color with others in the construction of II become 
uniquely colored in the construction of Il£. With a little extra care, one can 
obtain a coupling of n^ n ^' e (resp. II) and its reduction U.g' E ' r (t) (resp. U r s ), so that 
there are fewer blocks contained in IT*™)' 6 (resp. II) than in ~\j^' e ' r (resp. H r s ) at 
all times. 

Note that Il£ is also a S-coalescent, and that its driving measure is 



Then E S (A) = E S (A S ) = (1 - <5) 2 H(A). Let be defined as in Q2Q), but corre- 
sponding to n£, 



Since z i— > e z — 1 + z is an increasing function on [0, oo) we have ipg(q) < i/)(q)- 
In fact, (1271) states that ipg(q) = t/>((1 — S)q), q > 0, hence 



It is perhaps not a priori clear why Il£ is a simpler process. We will soon see 
that because its H-measure is concentrated on A 5 , the criterion of [21] for CDI is 
sharp, and under an additional condition, its asymptotic speed can be found in a 
way analogous to [1]. 

The second variation is as follows: given realizations of (|6]) and (J7J) as before, 
define the color- joining IF (j stands for "joining") of II to be the partition valued 
process where all the blocks with integral color are immediately merged together 
into one block. As for the 5-reduction, one can obtain a coupling of II^' 6 (resp. II) 



S 5 (c?x) := (1 - 5) E(d^)l {xeA s } 




(27) 
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and its color-joining nM^i(t) ( r esp. IP), so that there are fewer blocks contained 
in n( n ) ,£ ' J (t) (resp. IP) than in r±( n )' e (resp. II) at all times. 

The coalescent IP should be a A-coalescent, with its corresponding from 
( |T5l) given by 

^{9)= / J , Ul ^ J 5(dx). (28) 



A 



£~r^ 



The slick point is that the right-hand side in ( 128]) may be infinite. The existence 
of the integral in fl28|) is equivalent to (cf. the condition ((R|) in the next section) 



J A E;=1^ 



Indeed, we have 




(29) 



The upper bound is just an application of (JzT|) . For the lower bound, assume that 
q > 1, the argument is simpler otherwise. Note that e _z — 1 + z > z 2 /2 — z 3 /3 > 
z 2 /10 for z < 5/4, and that e~ z - 1 + z > z/10 for 2 > 5/4. Substituting 
z = q^2iLi x i we arrive at ( 1291) . As a consequence, the right-hand-side in ( 1281) is 
finite for one g G M+ if and only if it is finite for all q G 1R+. 



3 Main results 
3.1 Regular case 

In this subsection assume that (in addition to H({0} U A/) =0) the measure 5 
satisfies the regularity condition 

H(rfx)<oo. (R) 



The complementary setting is discussed in Section [3^21 

Denote by (N^(t),t > 0) the number of blocks process for the S-coalescent 
(r±(t),t > 0), and recall definition (126]) . Regularity (JR} implies 



/ 2 S ( rfx ) <0 °' 

JA\A 1 ~ a l^i=l x i 
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for any a G (0, 1). In particular, in the PPP construction an atom (t, x) satisfying 
Xi > a appears at a strictly positive random (exponential) time. Therefore, if 

S a ((ix) = H((ix)l{^ Xi <a}, the S-coalescent and the S a -coalescent have the same 

small time behavior. 

As already indicated in Section 12.31 the central result of this paper is 

Theorem 10. If both 2({0} U A f ) = and (EP hold, then 



where oo/oo = 1. In particular, under these assumptions, the candidate speed is 
finite if and only if the S-coalescent comes down from infinity, which happens if 
and only if it is the true speed of CDI. 

The proof is postponed until Section HI The importance of condition (jR|) will 
become evident in view of Lemma [18J that implies Proposition [T7] (see also (133]) - 
( |3^|) ). which is an essential ingredient in the martingale analysis of Section HI It is 
interesting to note that (JR]) arises independently in the context of the color-joining 
construction in Section 12.61 

Remark 11. Once given Theorem [T0| by straightforward copying of arguments 
from jl], one could obtain the convergence of N a /vs in the LP sense for p > 1, as 
well as the convergence of the total length of the genealogical tree in the regular 
setting. 

Lemma 12. Under ICR]) the color-joining IP is a k-coalescent corresponding to 
ipi from l[28\) . If IT comes down from infinity, then IP comes down from infinity 
at least as fast as U, meaning that IP has the speed of CDI vi(t) determined by 



This is an immediate consequence of elementary properties of the Poisson point 
process n from fl6]). 




1, almost surely, 




where vi(t) < f=(t), for any t > 0. 



Proof. Consider the process n( n ) = lim^o IT n )' e from the PPP construction of II. 
It suffices to show that n( n ) J is a A-coalescent (started from a configuration of n 
blocks) corresponding to 
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From the coupling of II and IP , where II has at least as many blocks as IP at 
any positive time, it is clear that if II comes down from infinity, then also does 
IP . Alternatively, the reader can verify analytically that 

(V%)-Ve(<?))'>0, g>0. 

Since IP is a A-coalescent, we know that vl is its speed of CDI. Then again due to 
the above coupling of II and IP we conclude that vl(t) < v~,(t), for any t > 0. □ 

Remark 13. Note that all the S-coalescents with 5 of the form (T4]) are regular, 
and more generally, if 5 is supported on any "finite" subsimplex {x : = 0, Vfc > 
n} of A, then the corresponding S-coalescent is regular. In particular, the 5- 
coalescents featuring in the selective sweep approximation of [14, 25J are regular. 



3.2 Non-regular case 

Assume H({0} U Af) = as in the previous subsection. The setting where 



l(dx) = oo (NR) 



is more complicated, and the small time asymptotics for such S-coalescents is only 
partially understood. 

Due to observations made in the previous section, (1NRI) is equivalent to the 
fact that the integral in fl2"%j) diverges. 

Lemma 14. Under liNR\) the color-joining IP is a trivial process containing one 
block at all positive times. 

Proof. As for Lemma [121 consider the prelimit coalescents 11^ and their color- 
jommgs nWj. it is easy to verify that (INRI) implies instantaneous coalescence of 
any two blocks of IT"^', almost surely. Indeed, the rate of coalescence for a pair 
of blocks is given by the integral in fINRp . □ 

Remark 15. The last lemma holds even if II does not come down from infinity. 

The following illuminating example was given in [21]. Suppose H has an atom 
of mass l/2 n at 

x := ( y x 1 , . . . , x 2 n_ 1 , 0, . . .), 

where = l/2 n , i — 1, . . . , 2 n — 1, and n G N. Then ips(q) = Q(q log(g)) so t>= is 
infinite, but the corresponding II comes down from infinity. Due to Theorem [10] 
we see that (TRl) cannot hold in this case. 
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It is useful to consider a generalization as follows: for a sequence / : PJ — > (0, 1), 
let E have atom of mass l/2 n at 



x' 



.71 




where again x" = 1/2", z = 1, . . . , [f(n)2 n \ , and where we assume that [f(n)2 n \ G 
{1, . . . , 2 n ~ x — 1}, n > 1, so that S(A/) = 0. It turns out that again ipn(q) — 
0(glog(g)) (in fact, this asymptotic behavior is uniform in the above choice of /), 
while the integral in ((R|) (or ( 1NRI) ) is asymptotic to 



Due to Theorem [10] we see that as soon as the above series converges, the cor- 
responding S-coalescent does not come down from infinity. However, its color- 
joining will in many cases come down from infinity, for example if f(n) = n~ 2 , 
then rf)S(q) = 9(g 3/2 ). 

Proposition 16. Suppose that $NR\) holds and that the the corresponding (stan- 
dard) S-coalescent IT has an infinite candidate speed (or equivalently, that (2B) 
fails). Then for any 5 G (0,1), its 5-reduction U r s does not come down from 
infinity. 

Proof. If II does not come down from infinity, then U r s does not either, due to the 
monotone coupling of II and Tig. 

Even if II comes down from infinity, we have that 



and, as already observed, that J°° l/tpg(q) dq = oo (for one and then all a G 
(0, oo)), so due to Proposition 33 of [21], il£ does not come down from infinity. □ 

The last result and Lemma [TH indicate the level of opacity of the non- regular 
setting. Indeed, a S-coalescent II that comes down from infinity, but has infinite 
candidate speed and satisfies (INRI) . can be formally "sandwiched" between its 
corresponding IF and Il£, where 5 > is very small, however the lower bound IP 
is trivial, and the upper bound Il£ does not come down from infinity, so one gains 
no pertinent information from the coupling. 

To end this discussion, let us mention another class of frustrating examples. 
Suppose that Hi is a probability measure on A satisfying both ([Rj) and ( [22]) and 



£/(")■ 



Tl 
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denote by v± the speed of CDI for the corresponding Hi-coalescent. Let S 2 be a 
probability measure on A satisfying flNRj) . Define 

so that 5 satisfies both (1221 and (INRI) . Due to easy coupling, the S-coalescent 
comes down from infinity, and moreover 

lim sup -r^— = 1 • 

The martingale technique however breaks in the non-regular setting, and we have 
no further information on the small time asymptotics of iV = . It seems reasonable 
to guess that iV = is asymptotic to v±/2 as t — > 0. 

Remark [22] discusses an approach that might be helpful in resolving the ques- 
tion of speed for S-coalescents that come down from infinity in the non-regular 
setting. 

4 The arguments 

The goal of this section is to prove Theorem [10] (that is, Theorem [1]) by adapting 
the technique from [I]. 

As already noted, the function ip defined in f[2Ul) is strictly increasing and 
convex. Furthermore, it is easy to check that v'(s) = —ip(v(s)) where v = v-~ is 
defined in (|26|) . so that both v and |i>'| are decreasing functions. 

Due to the observation preceding the statement of Theorem [TO], we can suppose 
without loss of generality that supp(H) C A 3 / 4 (recall notation ( EH]) ). As in [I], 
this will simplify certain technical estimates. 

To shorten notation, write iV instead of N~. Note that the function v is the 
unique solution of the following integral equation 

/"* ib(v(r)) 

hg(v(t)) - logMz)) + / y \ )J dr = 0, V0 < z < t, (30) 

Jz v[r) 

with the "initial condition" v(0+) = oo. If H can be identified with a probability 
measure A on [0, 1] as in (j3J), then ( 130]) is identical to the starting observation in 
the proof of Theorem [TO] for A-coalescents (cf. proof of [lj Theorem 1). 

Indeed, the rest of the argument is analogous to the one from jl] for A- 
coalescents, the general regular S-coalescent setting being only slightly more com- 
plicated. The few points of difference will be treated in detail, while the rest of 
the argument is only sketched. 
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4.1 Preliminary calculations 

Assume that the given S-coalescent has a finite number of blocks at some positive 
time z. Consider the process 



M(t) :=log(N(t))-log(N(z)) + 



ip(N(r)) 
N(r) 



dr. t > z. 



Let uq > 1 be fixed. Define 

T m : = inf{s > : N(s) < n }. 



(31) 



It turns out that, under the regularity hypothesis ([R]), M(t A r no ) is "almost" 
(up to a bounded drift correction) a local martingale, with respect to the natural 
filtration (J r t ,t > 0) generated by the underlying S-coalescent process. 

Proposition 17. There exists some deterministic n e N and C < oo such that 



E[dlog(N(s))\F s 



N(s) 



+ h(s) ds 



(32) 



where (h(s), s > z) is an T '-adapted process such that sup s6 [ z zATn ] \h(s)\ < C , and 

{s<T no } < C ds, almost surely. 
Both estimates are valid uniformly over z > 0. 

Restricting the analysis to n larger than n is a consequence of the following 
estimate, whose proof is given immediately after the proof of the proposition. 
Recall defined in ffTT]) . When taking probabilities or expectations with respect 
to the joint law of (Y £ , £ > 1), we include subscript x to indicate the dependence 
of the law on x. Define 

oo / oo \ * 

5(X) :=5> 2 + £: 



i=l 



i=l 



Lemma 18. There exists n G N and Co < oo such that for all n > n and all 

x G A 3 / 4 , we have 



lot 



and 



n — 



W - h 



{Y e (n) >0}J 



— \ogn) 



XXi nx£-l + (l- xt) 



n 



< C 5(x), 



n — 



£=1 



(n) 



1 {y/" ) >o}) 



-logn <C 5(x) 
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Proof, [of Proposition [77| Since ([Rj) holds, it suffices to show that for each s > 0, 
we have on {N(s) > no} 



E(d\og(N(s))\F 8 ) | ?KiV(s)) 



and 



ds ' N(s) 

E([d\og(N(s))] 2 \F s ) = O 



\h(s)\=0 



A 3 / 4 



Z-/i=l x » 



A 3 / 4 



Eoo 2 '— '(^ X ) ) dS, 
8=1 X i 



(33) 



(34) 



where O(-) can be taken uniformly in s. Note that the finite integrals above are 
in fact taken over A, since S is supported on A 3 / 4 . 

Recall the PPP construction of Section 12.41 and fix n > uq. On the event 
{N(s) = n}, an atom carrying value x e A arrives at rate l/(X]i=i x ?) S(dx) ds, 
and given its arrival, log iV(s) = logra jumps to log(n — Y^LiO^} — l, y („) >0 ,)). 
Therefore, 



£(dlog(jV(s))|.F s ) = / £ x 

JA 



log 



n 



H(dx) <is. 



Due to Lemma [18] and ( |23l) we can now derive ( |33l) . 

To bound the infinitesimal variance on the event {iV(s) = n}, use the second 
estimate in Lemma [181 together with the fact 



E([dlog(N(s))} 2 \F s 
ds 



< I E» 

'A 



( n J2e=i( Y e } 1 {y/"'>o} ■ 



log 2 



n 



i=l x 



E(dx) 



Finally, note that both ([33]) and (1341) are uniform upper bounds over s. 



□ 



Proof, [of Lemma [THj The argument is almost the same as that for [I] Lemma 
19 in the A-coalescent setting. Since the regularity "dichotomy" is a consequence 
of some more complicated expressions (arising in the calculations) in the current 
setting, most of the steps are included. Abbreviate 



Z (n) . 



>0}' 



n 



and note that Z^ 1 ' is stochastically bounded by a Binomial(n, xi) random 
variable. Let 

T = T n := log (1 - ZW) . 
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Split the computation according to the event 

A n = {ZW < 1/2}, 
whose complement has probability bounded by 

exp {-n (i log 1 + i log ^) } = 2V' 2 d - PT'\ 

uniformly in p := — an d n , due to a large deviation bound (for sums 

of i.i.d. Bernoulli random variables). On A c n we have |T| < logn, and on A n we 
apply a calculus fact, | log(l — y) + y\ < 2 (i- y ) — V 2 ^ V e IA 1/2], to obtain 

\E[T] + £ [Z<">Uj | < (logn)P(^) + E [{Z^fl An ] . 

Since < 1, we conclude 

\E[T] + E[Z^}\ < (logn + l)P{A c n ) + E[{Z^) 2 ]. 

Note that |£^[T] +£'[Z^] | is precisely the left-hand side of the first estimate stated 
in the lemma. Due to the estimate (52) in the proof of [I] Lemma 19 we have 

(log n)P{A c n ) < (logn)2> n/2 (l ~ p) n/2 < Cp 2 < CS(x), 

for some C < oo, all p G [0, 1/4], and all n large. 

Until this point the argument is identical to the one for A-coalescents. The 
new step is verifying that 

E[(Z^) 2 ] < S(x). (35) 
It is easy to check (see for example [3j Corollary 18) that 

E[{Yt ] - l {yf ) >0} ) 2 ] < Cn 2 (x e f, (36) 

for some constant C < oo. For two different indices k,£, use Cauchy-Schwartz 
inequality together with the above bound to get 

\E[{Yt n) - l {y w >0} )(^ (n) - V(«) >0} )]l < y/C*n*{x k nxtf = Cn 2 x k x e . (37) 

One obtains §5§ from (j36])-(j37|) after rewriting E[{Z^) 2 ] as 



V £ k i^k 

The second estimate is proved exactly as in [3]. □ 
Remark 19. The expectation of the product of — l, y (») >{) , and — 
l{y (n) >o} can ^ e com P u t e d explicitly, and one can verify that its absolute value 
has the order of magnitude n 2 XkXe as Xk and (or) xg tend to 0. 
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4.2 Proof of Theorem [TO 



Part I. Suppose that a given regular E-coalescent starts from n blocks, where 
n G N is large and finite. In other words, consider the prelimit process II^. 

Recall Remark [6j Define a family of deterministic functions (v n , n G N) as in 
(I18p. where ip = ip^, and note that v n satisfies t> n (0) = n and 



log{v n {t)) - log(n) + / t!^Ml dr o. Vf ;> o 



(3? 



/ w n (r) 

It is easy to see that the following is true. 

Lemma 20. We have v n (t) < v n+1 (t) and \im n ^. 00 v n (t) = Vs(t), for each t > 0. 

For each n > n Q (where n Q is the parameter from Proposition H7|) define the 
process 



N n (t A r n ) 
Mn (t) ;= log ; "V + , 



^(iV n (r)) ^(w n (r)) 



dr, t > 0, 



where h — h n is given in (|32|) . and r™ Q := inf{s > : N n (s) < n } in analogy to 

Due to Proposition [T71 and (138 p . we know that M n is a martingale (note that 
Af„(0) = 0), such that 

E[(M n (s) - M n (u)) 2 \T s ] < C(s - u), 

uniformly over n > no and u, s such that s > u > 0. Fix any a G (0, 1/2). Doob's 
L 2 -inequality therefore implies 



P(sup \M n (t)\>s a )=0(s 1 - 2a ), 
te[o,s] 

where O(-) term is uniform over n > n . Due to Proposition [T71 the term 



(39) 



(At? 



h(r)dr 



is of smaller order O(s), again uniformly in n > n . Hence we obtain from 
that 



P ( sup 

te[o,s] 



^Ar-) + 



^(A^ n (r)) ^(w n (r)) 
7V n (r) v n {r) 



dr 



> s L 



0{s 



l-2a\ 
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Due to Lemma [7] and [I] Lemma 10, the last estimate implies in turn 



sup 

te[o,s] 



log 



v n (tAT% ] 



> 2s c 



(40) 



Assume that the corresponding regular (standard) 5-coalescent II comes down 
from infinity. Since N n (t) /* N(t), for each t > 0, r™ Q /*■ r no , almost surely, and 
since 

P(r no > 0) = 1, (41) 

we obtain due to Lemma [201 that the candidate speed v(t) := lim n v n (t) is finite 
for each t > 0. 

Conversely, if this S-coalescent does not come down from infinity, then it must 
be v(t) := lim n v n (t) = oo. 

Part II. Suppose that the H-coalescent from part I comes down from infinity. 
It is tempting to let n — >■ oo in (j4"0~l) in order to obtain 



P [ sup 

t€[0,s] 



log 



N(t A T n 



v(t A r no ) 



> 2s a 



0(s 



l-2a\ 



(42) 



However, this step would not be rigorous without additional information on the 
family of events in (I40p . indexed by n > n . An alternative approach is discussed 
next. 

From part I we know that the corresponding candidate speed is finite. Using 
this fact, a variation of the previous argument yields f l4"2]) . Define a family of 
deterministic functions (v x , x £ M) by 

v x {t) = v(t + x), t > —x, 

and note that each v x satisfies an appropriate analogue of fl30|) on its entire domain, 
more precisely, v x (—x+) = oo and 



log(v x (t)) - log(v x (z)) + 



v x {r) 



dr = 0, V - x < z < t. 



(43) 



Due to PH. one can assume that z > r nn . For each x > — z define 



M ZiX {t) := log 



N(t A r, 



v x (t A r no ) 
where h is given in (132]) . 



v x lz) 



tAr„ 



*P(N(r)) _ ip(v x (r)) 
N(r) v x (r) 



h(r) 



dr, t > z, 
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It will be convenient to consider for each fixed oOa process M Zy x-, where 
X G T z such that P(X > —z) = 1. Note that such M Z)X is adapted to the 
filtration (J>, r > z). More precisely, let X z be the random variable defined by 

N(z) =v(X z + z) =v Xz (z). 

It is easy to see that X z + z is decreasing to as z decreases to 0, and that 
therefore the following is true. 

Lemma 21. We have lim^o^z = ; hence linx^o vx z (t) = v(t) for all t > 0, 
almost surely. 

Due to Proposition [T7I and ( 1431) . we know that M Z) x z is a martingale (note 
that M ZrX ,(0) = 0), such that 

E[(M z , Xz (s) - M Zi xXu)) 2 \F s ] < C(s - it) , 

uniformly over u, s such that s > u > z. As in part I, we obtain 



P(sup \M Zi xM\ > s a ) = 0(s 



l-2a\ 



te[z,s] 



where O(-) term is uniform over z > 0. Again, due to Proposition [T71 the term 

tf\T nQ 

h(r)dr 

is of smaller order O(s), uniformly in z. Hence 



P \ sup 

te[z,s 



VxAt A T n ) 



N{r) 



dr 



0(s 



l-2a\ 



As before, due to Lemma Hand [1] Lemma 10, the last estimate implies 

N(t A r m 



P ( sup 

te[z,s] 



log 



Ar„ ) 



> 2s c 



0(s 



l-2a\ 



and therefore for any z' < z 



P ( sup 

te[z,s] 



log 



N(t A r, 



> 2s a = CKs 1 -^) 



Let — >• and use Lemma [211 and then let z — > to obtain (|42|) . This together 
with (HTj) shows that in this setting the candidate speed is the true speed of CDI. 
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Remark 22. As already mentioned, the above argument works only under the 
assumption (jR|). However, regularity is only needed in linking E(dlogN(t)\J r t ) to 
—ip(N(t))/N(t), and in uniformly bounding the infinitesimal variance of logiV(i). 
For irregular H-coalescents that have an infinite candidate speed, but also come 
down from infinity, a relation of similar kind 

E(d\ogN(t)\F t ) = -^^ + h(t) 

might be possible, where h is still a uniformly bounded process, and where ipi is 
an increasing, convex function satisfying Lemma [7] and 

< oo, a > 0. 



It is natural to guess that v\ : M + — > K+, determined by f™r t )dq/ij>i(q) = t, is 
then the speed of CDI. 
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